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I. INTRODUCTION 


The initial research on inverted pendulums began in 
the 1950s, designed by control theory experts at the 
Massachusetts Institute of Technology (MIT) in the United 
States based on the principle of rocket launch boosters. 
The inverted pendulum system, as the model foundation 
for shipborne radar, rocket launch systems, and satellite 
attitude control, has been the focus of many researchers in 
the past few decades. The research on inverted pendulums 
will tend towards more complex and in-depth studies. 


Inverted pendulum systems can be divided into linear 
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Abstract— This study aims to conduct control research on an inverted 
pendulum system using the Simulink simulation platform. The inverted 
pendulum system is a classic nonlinear dynamic system with important 
theoretical and practical applications. Firstly, establish a mathematical 
model of the inverted pendulum system, including the dynamic equation of 
the pendulum rod and the sensor measurement model. Subsequently, the 
PID (proportional integral differential) controller design method based on 
the inverted pendulum system and the fuzzy PID controller design methods 
were verified through simulation experiments. The ultimate goal is for the 
designed fuzzy PID controller to effectively stabilize the inverted pendulum 
system in the vertical position and achieve fast tracking of the target 
position. Simulation and experimental results show that compared to 
traditional PID controllers, fuzzy PID controllers can quickly stabilize the 
pendulum in the target position and have good practicality, stability, speed, 
and accuracy. Future research can further explore the application of other 
advanced control strategies in inverted pendulum systems, as well as their 


potential applications in practical engineering. 


inverted pendulums, planar inverted pendulums, composite 
inverted pendulums, etc. according to their composition. 
According to their complexity, they can be divided into 
primary inverted pendulum systems, secondary inverted 
pendulum systems, tertiary inverted pendulum systems, 
and multi-level inverted pendulum systems. The first-level 
inverted pendulum system consists of a driving motor, a 
conveyor belt, a pendulum rod, a small car, and a test 
bench [1, 2, 3 ]. The first-level linear inverted pendulum 
system is driven by an electric motor and is an unstable, 


nonlinear, single-input, double-output, strongly coupled 
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system [4, 5, 6, 14, 15]. It controls both the angle of the 
pendulum and the position of the trolley to be stable, and 
the steady-state errors of the trolley position and pendulum 
angle must be controlled within a small range. 

The control of the inverted pendulum is a difficult 
point in the study of inverted pendulum control, and there 
have been many studies on the inverted pendulum, which 
are basically based on the assumption that the trolley track 
of the inverted pendulum system is sufficiently long [7]. 
With the development of technology, new control methods 
are constantly emerging, and people use inverted 
pendulums to test whether new control methods can 
handle multivariable, nonlinear, and absolute instability. 
The inverted pendulum has become an ideal experimental 
method for testing the effectiveness of control strategies 
[8]. This article focuses on a first-order inverted pendulum 
system. Firstly, a mathematical model is established using 
the knowledge of Newtonian mechanics, and then a 


simulation model of the inverted pendulum system is 


Data 
acquisition 
card 


International Journal of Advanced Engineering Research and Science, 10(12)-2023 


established using the Simulink module of MTALAB [9, 10, 
11, 13]. Finally, by comparing and analyzing the curves 
and parameters of the established traditional PID controller 
and fuzzy PID controller, this study is trying to find out 
whether the fuzzy PID control method is better or not than 
the ordinary PID control method in terms of stability and 


speed. 


II. ESTABLISHING MATHEMATICAL MODEL 

The working principle of a first-order linear inverted 
pendulum (Figure 1) is that when the data acquisition card 
transmits the collected data from the rotary encoder to the 
computer and compares it with the set value. The deviation 
is processed through some calculation, and a control law is 
issued to control the motor to make the pendulum swing 
left and right into the stable range, thereby achieving the 
pendulum to stand upright and not fall, as well as 


self-swing [12]. 


Inverted 
pendulum 


Fig.1 Working Principle Diagram of a First Level Linear Inverted Pendulum 


Because establishing a mathematical model of the 
system is the foundation for studying control methods, the 
first step is to model the inverted pendulum system in this 
paper. And mathematical modeling is carried out using the 
Newtonian mechanics method to obtain the state-space 
equation of the system and prepare for the subsequent 
controller design and simulation. 

The model parameters of the inverted pendulum 
system are the pendulum mass m\=0.109kg, the trolley 
mass m2=1.096kg g, the angle between the pendulum and 
the vertical direction 0 (rad), the distance from the center 
of the swing rod to the car /=0.25m, the distance the car 


moves x (m), the force applied to the car f (N), the friction 
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coefficient of the car rf =0.1N/m/sec, the inertia of the 
swing rod 1=0.0034kg*m2, and the gravitational 
acceleration g=9.8N/ m2. The physical model diagram of 
the inverted pendulum system is shown in Figure 2. A 
detailed decomposition of various forces acting on the 
pendulum and trolley using Newtonian mechanics 
methods is shown in Figure 3. P and N are set as the 
components of the interaction force in the vertical and 
horizontal directions during the movement or stability of 


the car. 
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Fig.2 Physical Model of Inverted Pendulum System 


N 


Fig.3 Force Analysis Diagram of the Small Car and 
Swing Rod 


According to Figure 3, first analyze the force in the 
horizontal direction. 


m: = F — rž — N (1) 
2 . 
N = m, <5 (x + Isin 0) = m, # +m,16 cos 0 — 


m,102 sin 0 (2) 

Substituting equation (2) into equation (1) yields the 
equation of motion in the horizontal direction 

(mı + m2) + r¢X + m,16 cos @ — m,16? sin@ = 
F (3) 


In the vertical direction, as the car is used as the horizontal 


plane, only the force acting on the swing rod needs to be 


analyzed to obtain 
a? sa 
P—mg= m Fall cos 0) = —m,16 sin 8 — 


m,l6 sin ð —m,l02cos@ (4) 


The torque balance equation is 


—Pl sin 0 — Nl cos 0 = IÖ (5) 


The second equation of motion can be obtained by using 


equations (2) and (4) 


(I+ ml*)¢ -mgl = m, lä (6) 
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Set angle 0= n+ ọ, and ọ after being converted to radians, 


it is much less than 1 rad, that is p <« 1. So it can be 


2 
simplified as cos@ = —1, and sin@ = —9, (a) E 


Using u instead of the input force F, equations (3) and (6) 
can be simplified as: 
l (+m l? )ö -m glo'= m lë 
(m, + m)X +r; x -—m,lg@ =u 
(7) 
Perform a Laplace transform on the above equation, which 
will be 
(I + m,17)®(s)s* — m,gl@(s) = m,1X(s)s? 
(M + m)X(s)s* + 7X(s)s — m, l®(s)s* = U(S) 


(8) 
The first equation in the above equation can be written as 
follows: 

(S) _ mls? 


(9) 


X(S) (+m ,12)—m gl 
Substitute equation (9) into the equation with the control 
input in equation (8), which will be 


OR 
uS) 


maıl 
1°62 
w 


st+ 


"fd+m4 12) 3 (mı+m2)Mmıgl 2 rpms | 
w w w 


(10) 
wherein, w = [(m, + m3) (I + mI”) — (m,1)?]. 
Using equation (7) for X and@ solve, which will be 


x=x 


-(1+m41° )rf mı? gl? (1+m11?) 


z? t 


(mı+mz2)I +m m31? (mi+mz)l +m mz 


o=9 


mylr¢ m,gl(m,+mz) mıl 


OEE S 
(mi+mz2)I+mım31? 


(mı +m2)I +m m31? (mı +m2)I +m m31? 


(1) 


Compiled into a standard state space equation as 


1 0 
-(1+m41° )rf m4? gl? 


(my+mMg)l+mymz2l2 (m1+mz2)I +m; m31? 


SS RR 


x 
x 
p 
7) mylr ¢ mıgl(mı+m2) 


O h O O 


0 
0 
10 


(mı +m2)I +m m31? 


(1+m11?) 
(mi+m3)I+mım3l? u (12) 


maıl 


(mi+mz)I+mm3l? 
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x 
y=(J=[! 0 0 oj|*ž 
0 0 0 1 œj? 
ọ 
If the pendulum is regarded as a rod with uniform mass, 
then the inertia of the pendulum is 


I= im? (13) 


Substituting equation (13) into the first equation of 
equation (7) yields 


(Emi? + ml?) ë -m glo = milk (14) 

Simplify the above equation, which will be 

= p+7ä (15) 

Let the system state space equation be 

y = AX + Bu (16) 
Y = CX + DU 


Let X=[x x 6 6],u =X, The following state 
space expression can be obtained 


x] 0 1 0 O7;x7 pO 
äl Jo 0 o offf Jij, 
g|—~|9 0 0 1]}Jp}t}ol¥ (17) 
əl lo o o {le} |= 
P 4l 4l 
x 
_ x7 71 0 0 Q|ž 
Y=lal=lo o 1 allo 
? 


IHI. CONTROLLER DESIGN AND FUZZY LOGIC 
ESTABLISHMENT 

Firstly, traditional PID is used to design PID 
controllers for the output displacement and output angle of 
the system. Through parameter tuning and control, the 
spatial-state equation output reaches a stable state. When 
studying the displacement of a small car, the spatial state 
control equation of the system is obtained by inputting the 


car parameters into equation (17) as follows: 


010 0 0 
_|0 0 0 o0 i lee 
seas | o 0o 1 a ee 

0 0 0 29.4 0.75 (18) 


Y=[xJ=[1 0 0 OlX 
We first use the traditional PID control method to 


control the spatial-state equation of displacement output 
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and simulate it using Simulink in Matlab. We then adjust a 
suitable set of PID parameters through parameter tuning 
principles. The input is the unit step signal, and the PID 
parameter is K,= 70, K:= 0.7, Kp= 6 the system reaches 
steady state at time T=0.08S. 

The fuzzy controller designed here is a two-input, 
three-output fuzzy PID controller, and the two output 
signals jointly use a fuzzy rule. The input is the error e= 
r(k)-y(k) and error change rate ec = e(k)-e(k-1) between 
the given value and the actual value of the car 
displacement or swing rod angle, and the output is the 
corrected values 1K,, AK;, and4Kp of the PID parameters. 

In this design, the basic domain of error e is taken as 
[-5, 5], and the basic domain of error change rate ec is 
taken as [-5, 5]. The domain of the output variablesAK,, 
AKi, andAKp are taken as [-3, 3]. 

In order to obtain the input of the fuzzy controller, it 
is necessary to fuzzify the precise quantity, that is, 
multiply the input quantity by the corresponding 
quantization factor, and convert it from the basic domain 
to the corresponding fuzzy domain. The quantization 
factor of error e isa-=0.8, and the factor of error change 
rate ec iS Oec=0.2. The control quantity obtained through 
the fuzzy control algorithm is a fuzzy quantity that needs 
to be multiplied by a proportional factor and converted 
into the basic domain. When the output variable is the 
displacement of the car, the scaling factor of4K,,4K,4Kp 
areaakp =Gak1=1,aakp =-5. When taking the output variable 


swing angle, the scaling factor of AK,4K,4Kp 


areaakp=200,aaxi—1,aaxp=30 

Divide the fuzzy domain of input variables(e, ec) and 
output variables (4K,,4K;,A4Kp) into 7 fuzzy subsets, 
namely NB, NM, NS, ZO, PS, PM, PB representing 
negative big, negative medium, negative small, zero, 
positive small, positive medium, and positive big, 
respectively. The membership functions of input variables 
and output variables both adopt triangular membership 


functions, as shown in Figure 4. 
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NB NM NS 0 PS PM PB 


Fig.4 The Membership Function of e, ec, AK AK, AKp 


The fuzzy rule statement used by the output variable fuzzy controller is as follows: 
“if E isaand EC isfthen U isy” Wherein,a,f,yboth represent the fuzzy sets corresponding to each variable. 
Based on the impact of PID parameters on system performance, parameter tuning principles, expert experience, and 
cognition, 49 control rules were obtained after processing, as shown in Tables 1-3. 
Table 1 Fuzzy Rule Table ofAK, 


EC 
NB PB PB PB PB PM PS 
NM PB PB PB PB PM O 
NS PM PM PM PM O PS PS 
O PM PM PS O NS NS NM 
PS PS PS O NS NM NM NM 
PM PS NS NM NM NM NB 
PB O NM NM NM NB NB 


i NB NM NS O PS PM PB 
EC 

NB NB NB NM NM NS 
NM NB NB NM NS NS 

NS NB NM NS NS O PS PS 

O NM NM NS O PS PM PM 

PS NM NS O PS PS PM PB 

PM PS NM PM PB PB 

PB PS PM PM PB PB 
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Table 3 Fuzzy Rule Table ofAKp 


E 
z NB NM NS O PS PM PB 
EC 
NB PS NS NB NB NB NM PS 
NM PS NS NB NM NM NS 
NS O NS NM NM NM NS 
O O NS NS NS NS NS 
PS O O O O O O 
PM PB PS PS PS PS PS PB 
PB PB PM PM PM PS PS PB 


This design system uses the Mamdani inference 
method to perform fuzzy inference on the established 
fuzzy rules in order to obtain control variables. Meanwhile, 
using the center of gravity method to solve the fuzziness of 
language expression, thus obtaining the exact value 
of4K,,4K,AKp. In addition, the values obtained through 
fuzzy reasoning and deblurring are multiplied by the 
corresponding scaling factors to obtain the incremental 
adjustment values of PID parameters, which are then 
substituted into equations (19) - (21) to obtain the control 


parameters of the PID controller. 


Kp = Kpo + A Kp (19) 
K, = Ko + AK; (20) 
Kp = Koy + A Kp (21) 


Before establishing a fuzzy logic controller, the Tool 
Box parameters of the fuzzy logic system need to be set in 
Simulink (Figure 5). Then use the membership function 
parameter setting process of Tool Box (Figure 6) to 
establish the membership function of e, ec,4K,,4K,4Kp. 


“& Fuzzy Logic Designer: Fuzzy_PID 
File Edit View 


= oO x< 


Fuzzy_PID 


(mamdani) 


Fuzzy_PID 
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System "Fuzzy_PID": 2 inputs, 3 outputs, and 49 rules 


Fig.5 Fuzzy Logic System Tool Box 
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| “& Membership Function Editor: Fuzzy_PID 


File Edit View 


FIS Variables 


Membership function plots 
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Fig.6 Membership Function of Fuzzy Logic System 


After setting the above parameters, they can be 
added to the rule editor, and all output surfaces of the 
fuzzy inference system can be observed in the output 


surface observer. 


IV. COMPARISON BETWEEN TRADITIONAL 
PID CONTROL AND FUZZY PID CONTROL 
Through Simulink simulation (Figure 7 and Figure 


8), it was found that, under appropriate parameters, 


Gain11 


Gain12 


Derivative2 


traditional PID controllers have good control effects on the 
control object that can assume an accurate mathematical 
model. They can meet the requirements of system control 
accuracy and rise time, but there are problems such as long 
adjustment times. The fuzzy PID controller can adjust the 
PID parameters of the system in real time. By doing so, the 
PID parameters can be more suitable for the control 
requirements of the system, resulting in better control 
effects than traditional PID controllers (Figure 9). 


x = Ax+Bu 
y = Cx + Du 


State-Space4 


Add4 
Mux1 Scope7 


Fuzzy pid2 


Fig.7 Simulink Simulations of Traditional PID and Fuzzy PID Control Systems (Swing Angle) 
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Fuzzy Logic 


Derivative Controfert 


Producti 


Product2 


Simulation Comparison of 
Traditional PID Control and Fuzzy PID Control 


0.5 PID Controller 
‘ Fuzzy PID Controller 
0 0.1 0.2 0.3 0.4 


0.5 


Times(s) 


0.6 0.7 0.8 0.9 1 


Fig.9 Response Curves of Traditional PID and Fuzzy PID Controller Systems (Swing Angle) 


V. CONCLUSIONS 

The inherent first-order linear inverted pendulum 
system is a typical single-input, multiple-output nonlinear 
system. This article analyzes the motion of the inherent 
first-order linear inverted pendulum system, establishes a 
mathematical model of the first-order linear inverted 
pendulum system using the Newtonian mechanics method, 
and designs PID controllers and fuzzy PID controllers to 
control the first-order linear inverted pendulum system 
separately. 

By adjusting their proportional constant and integral 
constant, the three parameters of the differential constant 
enable the first-order linear inverted pendulum system to 
quickly and accurately reach a stable state, ultimately 
achieving a stable system. The superiority of fuzzy PID 
controllers over PID control was verified through Simulink 
simulation. Under the action of a fuzzy PID controller, the 
system has fast response speed, short adjustment time, and 


high steady-state accuracy, achieving the expected goals. 
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